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Abstract 

This paper presents a polynomial-time 1/2-approximation algorithm for maximizing nonneg¬ 
ative fc-submodular functions. This improves upon the previous max{ 1 /3,1 /(l+a)}-approximation 
by Ward and Zivny [15], where a = max{l, sj (fc — 1)/4}. We also show that for monotone fc- 
submodular functions there is a polynomial-time fc/(2fc — l)-approximation algorithm while for 
any e > 0 a ((fc + l)/2fc + ^-approximation algorithm for maximizing monotone fc-submodular 
functions would require exponentially many queries. In particular, our hardness result implies 
that our algorithms are asymptotically tight. 

We also extend the approach to provide constant factor approximation algorithms for maxi¬ 
mizing skew-bisubmodular functions, which were recently introduced as generalizations of bisub- 
modular functions. 


1 Introduction 

Let 2 V denote the family of all the subsets of V. A function g : 2 V —>• M is called submodular if it 
satisfies 

g(z i) + g(Z 2 ) > g{Z\ u Z 2 ) + g{Z\ n Z 2 ) 

for every pair of Z\ and Z 2 in 2 V . Submodular function maximization contains important NP- 
hard optimization problems such as max cut and certain facility location problems. It is known 
to be intractable in the standard value oracle model, and approximation algorithms have been 
studied extensively. In particular, Feige, Mirrokni, and Vondrak [6] have developed constant factor 
approximation algorithms for the unconstrained maximization of nonnegative submodular functions 
and shown that no approximation algorithm can achieve the ratio better than 1/2. Buchbinder, 
Feldman, Naor, and Schwartz [2] provided much simpler algorithms that substantially improve the 
approximation factor. In particular their randomized version, called the randomized double-greedy 
algorithm, achieves the factor of 1/2, which is the best possible in the oracle value model. 
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In this paper we shall consider the maximization of nonnegative fc-submodular functions which 
generalizes the submodular function maximization. Let (k + 1)' := {(Xi, • • • ,X k ) \ Xi C V (i = 
1,..., k),Xi fl Xj = 0 (i j)}. A function / : (k + l) v —» K is called k-submodular if, for any 
x = (Ah,..., X k ) and y = {Y\ ,..., Y k ) in [k + l) v , we have 


f(x) + f{y) > /(* u y) + f{x n y ) 

where 


x\~\ y 
xUy 

X-submodular functions were first introduced by Huber and Kolmogorov [9] as a generalization 
of bisubmodular functions , which correspond to 2-submodular functions in the above notation. 
Examples of bisubmodular functions include the rank functions of delta-matroids and the cut 
capacity functions of bi-directed networks, and the minimization problem has been extensively 
studied [8, 12]. Examples of /c-submodular functions will be explained later. 

Ward and Zivny [15] and the present authors [11] independently observed that algorithms for 
submodular function maximization due to Buchbinder, Feldman, Naor, and Schwartz [2] can be 
naturally extended to bisubmodular function maximization. In particular the randomized double 
greedy algorithm for submodular functions can be seen as a randomized greedy algorithm in the 
bisubmodular setting and it achieves the best approximation ratio 1/2. Ward and Zivny [15] 
further analyzed the randomized greedy algorithm for fc-submodular function maximization and 
proved that its approximation ratio is 1/(1 + a), where a = rnaxjl, y/(k — 1)/4}. They also gave a 
deterministic 1/3-approximation algorithm. 

In this paper we shall present an improved 1/2-approximation algorithm for maximizing k- 
submodular functions. Our algorithm follows the randomized greedy framework as in [11, 15], and 
the main idea is the use of a different probability distribution derived from a geometric sequence 
at each step. 

By extending the argument by Feige, Mirrokni, and Vondrak [6] we also show that for any 
e > 0, a ((k + l)/2fc + ^-approximation for the X-submodular function maximization problem 
would require exponentially many queries, implying the tightness of our result for large k. In fact, 
our inapproximability result holds for a much restricted class of monotone fc-submodular functions, 
where a /c-submodular function is said to be monotone if f(x) < f(y) for any x = (X\,... ,X k ) 
and y = (Yf,..., Y k ) in (k + 1)' with C F, for 1 < i < k. On the other hand we show that 
there is a k/(2k — l)-approximation for monotone fc-submodular functions. In particular it attains 
an approximation ratio of 2/3 for bisubmodular functions. 

In order to understand the relation between fc-submodular function maximization and other 
maximization problems, it is useful to understand characteristic properties of /c-submodular func¬ 
tions, called orthant submodularity and pairwise monotonicity. To see them, define a partial order 
Y on {k + l) v such that, for x = (Xi ,..., X k ) and y = (Yi,..., Y k ) in {k + l ) v , x A y if Xi C Yj 
for every i with 1 < i < k. Also, define 


:= (X 1 r\Y 1 ,...,X k nY k ) 


X! U Yi \ (J Xi U Yi ., X k U Y fc \ (J Xi u Yi 


v^/1 




A e,if(x) = f(Xi, ..., Xi-t,Xi U {e}, X i+1 ,. 


..,X k )-f(X 1 ,...,,X k ) 
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for x € (k + l ) v , e ^ (J^ =1 Xj, and i € [A:], which is a marginal gain when adding e to the z-th 
component of x. Then it is easy to see that the fc-submodularity implies the orthant submodularity: 


A e ,if(x) > A e ,if(y) (x , y € (k + l) y with x<y,e$ U Yj, and i € [A]), 

ie[fc] 

and the pairwise monotonicity: 

A e ,if(x) + A e jf(x) >0 (* € (fe + l) y , e ^ |J 2Q, and i, j € [fe] with z / j). 

£e[fc] 


Ward and Zivny [15] showed that these properties indeed characterize A-submodular functions, 
extending the corresponding result for bisubmodular functions [1]. 

Theorem 1.1 (Ward and Zivny [15]). A function f : (k + 1)' —>• R is k-submodular if and only if 
f is orthant submodular and pairwise monotone. 

The fc-submodular function maximization problem is closely related to the submodular function 
maximization with a partition matroid constraint. Consider a partition {U\, ..., U n } of a finite set 
U such that \U t \ = k and a partition matroid on U such that I C U is independent if and only if 
1 1 DU,] < 1 for every 1 < i < n. By identifying each Ui with [fc], one can identify each independent 
set I with an element x of (k + l) 1- , where V = {1, ... ,n}. Therefore, for a given submodular 
function g : 2 U —>• M + , its restriction to the family of independent sets can be considered as a 
function from (k + l) v to M + satisfying orthant submodularity. In general, if g is monotone, the 
submodular function maximization with a matroid constraint admits (1 — l/e)-approximation [4], 
which is known to be best possible in the value oracle model [13]. On the other hand, when g is 
non-monotone, the current best approximation ratio is 1/e [7] for general matroids, and deriving 
the tight bound is recognized as a challenging problem even for uniform matroids (see [3]). The k- 
submodular function maximization is in between: it admits 1/2-approximation whereas it assumes 
pairwise monotonicity, which is strictly weaker than monotonicity. 

It is also worth mentioning that in the A-submodular function maximization there always exists 
a maximizer which is a partition of V (c.f. Proposition 2.1), which corresponds with a base in the 
partition matroid. Vondrak [14] showed that, under a matroid base constraint, any (1 — \jv + e)- 
approximation requires exponentially many queries for any e > 0, where v denotes the fractional 
packing number (see [14] for the definition). One can easily show that v = k in our case, and hence 
this general result does not give a nontrivial bound for large k. 

We should also remark that, in the /c-submodular function maximization problem, function 
values are determined over (k + I) 17 and hence over the independent sets in the corresponding 
submodular function maximization with a partition matroid constraint. It is not in general true that 
such a non-negative (monotone) function can be extended to a non-negative (monotone) submodular 
function over 2 U . 

An important special case of the submodular function maximization with a partition matroid 
constraint is the submodular welfare problem. In the submodular welfare problem, given a finite 
set V and monotone submodular functions gt : 2 V —>• for 1 < i < k, we are asked to find a 

partition {X\,... ,X^\ of V that maximizes Yli=i 9i{Xi)- Feldman, Naor and Schwartz [7] gave a 
(1 — (1 — 1/A) fc ^approximation approximation algorithm, which is known to be best possible in the 
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value oracle model [14]. Now, consider h : (k + 1)^ —>• M+ given by 

k 

h(X 1 , ...,X k ) = Y^ 9i(Xi ) ((* 1 , ...,**)€(* + 1) V )- 

Then the submodularity and the monotonicity of g t imply the orthant submodularity and the 
pairwise monotonicity of h, and hence h is monotone /c-submodular by Theorem 1.1. Thus the 
monotone /c-submodular function maximization generalizes the submodular welfare problem. In fact 
we will show that the approximation algorithm by Dobzinski and Schapira [5] for the submodular 
welfare problem can be extended to the monotone case. 

A similar construction gives another interesting application of the /c-submodular function max¬ 
imization. For a submodular function g : 2 1 —>• M_j_, define b! : (k + \) v —» M + by 

k 

h'(x 1 ,...,x k ) = J2g(x i ) ((x 1 ,...,x k )e(k + i) y ). 

2=1 

The resulting h! satisfies orthant submodularity but may not satisfy pairwise monotonicity in 
general. However if g is symmetric (i.e., g(X) = g(V \ X) for X C V) it turns out that h! 
is pairwise monotone and thus it is /c-submodular by Theorem 1.1. Therefore, for a symmetric 
submodular function g, our algorithm gives a ^-approximation for the problem of finding a partition 
{Xi ,..., X k } of V that maximizes X^i=i /7(A) • Note that this problem generalizes the Max k- cut 
problem. 

As another extension of the bisubmodularity, Huber, Krokhin, and Powell [10] have introduced 
the concept of skew-bisubmodularity. For a € [0,1], a function / : 3^ —> R is called a-bisubmodular 
if, for any x = (X\, X 2 ) and y = (Yi, Y 2 ) in 3 V , 

f(x) + f(y ) > f(x n y) + af(x Uy) + (1- a)f(xUy), 


where 

xUy = (Ai U Yi, (X 2 U Y 2 ) \ (Xi U Hi)). 

A function / : 3 V —>• R is called skew-bisubmodular if it is a-bisubmodular for some a € [0,1]. 

We show that a randomized greedy algorithm provides an approximate solution within the 
factor of f° r maximizing an a-bisubmodular function. This means that the double greedy 

algorithm of Buchbinder et al. [2] relies on a symmetry of submodular functions. Combining this 
with another simple algorithm, we obtain an approximate algorithm whose approximate ratio is at 
least for any a € [0,1]. This result has been included in our previous technical report [11], but 
not in a reviewed article. 

The rest of this paper is organized as follows. In Section 2, we present our approximation algo¬ 
rithms for the /c-submodular function maximization. In Section 3, we discuss the inapproximability. 
In Section 4 we analyze a randomized greedy algorithm for maximizing a-bisubmodular functions, 
and then we present an improvement that leads to a constant-factor approximation algorithm. 

2 Approximation algorithms for /c-submodular functions 

In this section we give approximation algorithms for the /c-submodular function maximization prob¬ 
lem. To analyze /c-submodular functions it is often convenient to identify (/c + 1)^ as {0,1... , k} v , 
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that is, the set of |l/|-dimensional vectors with entries in {0,1,..., k}. Namely, we associate 
(X,..., Xk) € (k + 1)' with x € {0,1,... , k} v by Xi = {e € V \ x(e) = i} for 1 < i < k. 
Hence we sometimes abuse notation, and simply write x = (X,... ,Xk) by regarding a vector x 
as a subpartition of V. 

For x € {0,1,..., k} v , let supp(a;) = {e € V \ x(e) ^ 0}, and let 0 be the zero vector in 


{0,1,... , k} v . 


2.1 Framework 

Our approximation algorithms are obtained from the following meta-framework (Algorithm 1) for 
maximizing /c-submodular functions by changing the probability distributions used in the frame¬ 
work. 

Algorithm 1 

Input: A nonnegative /c-submodular function / : {0,1,... , k} v —> M+. 

Output: A vector s. 

S i — 0. 

for each e € V do 

Set a probability distribution p over {1,..., k}. 

Let s(e) € {1,... , k} be chosen randomly, with Pr[s(e) = i] = pt for all * € {1,... , k}. 

return s 


The approximation algorithms for bisub modular functions [11] and more generally for /c-submodular 
functions [15] are specializations of Algorithm 1, where the probability distribution is chosen to be 
proportional to its marginal gain. 

We now evaluate the quality of the solution of Algorithm 1 by applying the analysis in [11, 15]. 

We first remark the following key fact (see [11, 15] for the proof). 

Proposition 2.1. For any k-submodular function / : (/c + 1)' —>• there exists a partition ofV 

that attains the maximum value of f. 

We also need the following notation, which will be used throughout this section. Let n = \V\. 

By Proposition 2.1 there is an optimal solution o with supp(o) = V. Let s be the output of 
the algorithm. We consider the j-th iteration of the algorithm, and let be the element of V 
considered in the j-th iteration, pf be the probability that i- th coordinate is chosen in the j-th 
iteration, and be the solution after the i-th iteration, where s = 0. Also for 0 < j < n 
let cf 'H = (o U s^) U s^\ that is, the element in {0,1,..., Zc} v obtained from o by replacing the 
coordinates on supp(s^) with those of s^\ and for 1 < j < n let 1 -* = (oUs^)Us^ _1 ( that is, 

the one obtained from by changing o^(e^) with 0. Also for i € [k] let yf' > = A e p) 
and let = A e yj Due to the pairwise monotonicity, we have 



( 1 ) 

(2) 


Also from s^ A t^\ the orthant submodularity implies 



( 3 ) 
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Applying the analysis in [11, 15], we have the following. 

Lemma 2.2. Let c G M+. Conditioning on suppose that 

k k 

K-* } - a i } )p { i ] < c E i/P V? } ) (4) 

i=l i=l 

holds for each j with 1 < j < n, where i* = o(e^). Then E[/(s)] > 

Proof. Conditioning on s^- 1 ), we have E[/(o^ -1 ^) — f(o^)] = J2i( a P ~ a P)pP an< 4 E [/( s< ^) — 

f(s^~ 1 ' ) )} = YhivP pP■ Hence, by (4), we have E[/(oh’ -1 )) — /(o^)] < cE[f(s^) — /(s^ _ H)] 
(without conditioning on s ^~ 1 )). Note also that c/ 0) = o and o ( - n -* = s by definition. Hence 

n 

f(o) - E[/(s)] = E [/(o°' _1) ) - /(° 0) )] 

j = l 

n 

< c (]Te [/(^) -/(^'- 1) )D 

i=i 

= c(E[/(a)] - /(0)) < cE[/(s)], 


and we get the statement. □ 

2.2 A ^-approximation algorithm for non-monotone fc-submodular functions 

In this section, we show a polynomial-time randomized ^-approximation algorithm for maximizing 
fc-submodular functions. Our algorithm is described in Algorithm 2. 

Theorem 2.3. Let o be a maximizer of a k-submodular function f and let s be the output of 
Algorithm 2. Then E[/ (s)] > \ f (o). 

Proof. By Lemma 2.2 it suffices to prove (4) for every 1 < j < n for c = 1. For simplicity of the 
description we shall omit the superscript (j) if it is clear from the context. Our goal is to show 

^2 (Vi + a i)Pi > a i*, ( 5 ) 

1 <i<k 

which is equivalent to (4) with c = 1. Recall that y ? ; + yy > 0 and ai + a*/ > 0 for i, i! G [k] with 
i / if, and y t > ai for i G [k] (c.f. (1), (2) and (3)). 

If i + < 1, then we need to show a\ + y\> ai*. Since yi + yv > 0 for i, if G [k], we have y\ > 0. 
Hence a\ + y± > ai* holds if i* = 1. If i* ^ 1 then 0 > yi* > ai* , and hence a\ > 0 by a± + a** > 0. 
This implies ai + yi > 0 > a^. 

If i + = 2 , we need to show (or + yi)yi + (a 2 + 272)2/2 > a**(yi + y 2 ). Now (ai + y\)y\ + (a 2 + 
2/2)2/2 = aiyi + a 2 y 2 + (yi - V2 ) 2 + 2 yiy 2 > aiyi + a 2 y 2 + 2 yiy 2 . If i* = 1 , then aiyi + o 2 y 2 + 
2 yiy 2 > ai(yi + y 2 ) + (a 2 + ai)y 2 > ai(yi + y 2 ) as required. By a symmetric calculation the 
claim follows if i* = 2 . If i* > 3, then 0 > yi* > a**, and hence ai > 0, a 2 > 0. We thus have 
(ai + yi)yi + («2 + y2)V2 > o > a*.(yi + y 2 ). 
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Algorithm 2 

Input: A nonnegative fc-submodular function / : {0, !,■■■, L} 1 —>• M+. 
Output: A vector s G {0,1, ... , k} v . 

Si— 0. 

for each e € V do 


0 + 


if 


Hi A e> j/(s) for 1 < i < k. 

Assume yi > y 2 > ■ ■ ■ > yk- 

{ the maximum integer i such that yi > 0 
0 

< 1 then 

_fl if * = 1 

1 0 otherwise 

else if = 2 then 

“^< 1 , 2 } 

I 0 otherwise 

else 


Pi 


(1 < i < k). 


(1 <i< k ). 


if y i > 0, 
otherwise. 


Pi <~ 



i 


if i < i + — 1 
if i = i + 


(1 < i < k). 


0 otherwise 

Let s(e) € {1,... , k} be chosen randomly, with Pr[s(e) 

return s 


i] = pi for all * G {1,, k}. 


Hence assume i + > 3. Note that 


Vi > yi* > ai* for i < i*. (6) 

Let r G argminjaj | i G [A:]}. Such r is unique if dr ^ 0* 

If r = i* , we have Yli a iPi — a i*^2,iPi) = a i* ■ Since YliViPi — 0) (5) follows. Hence we assume 
r ^ i*. 

If i* > i + , we have E; y t p t = Ei<i+ ViPi > E,:<i+ a i*Pi = a i* by ( 6 ) and E; a iPi = Ei^r a iPi + 
a r p r > 0 by Yli^rPi — Pr anb dj + a r > 0 for i / r. Therefore (5) holds. We thus assume i* <i + . 
Now we have 


22(Vi + ®i)Pi > 22 ai * Pi + 22 aiPi + 22 aiPi 


i<i* 


i>i* 


y, CLi*pi + 2 a,i*pi* ) + j y aiPi + y aiPi + a r p r 


Ki<i 


For the first term we have 


^ ^ di*Pi + 2 cii*pi* 




(7) 
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Hence it suffices to show that the second term of (7) is nonnegative. This is trivial if a r > 0. Hence 
assume a r < 0. Since i* < i + , we have 

X pi + X pi = Q) + X pi= pi *+ X pi = 1 - p r- ( § ) 

i>i* 

Therefore, if r < i*, we get 

^ ^ 0 >iPi T ^ ) &iPi T a r p r ^ I p r ^ ) Pi ^ ) Pi I = ®r(Pr (1 Pr )) = a r (2p r 1) -L 0, 

i>i* i^r,i* \ *>** ) 

where the first inequality follows from Oj + a r > 0 for i / r, the second equality follows from ( 8 ), 
and the fourth follows from a r < 0 and p r < 1/2. Hence we further assume r > i* . Then p r < 1/4 
by r / 1 and i + > 3. Hence, by a r < 0, 


X aiPi + X aiPi + arPr - X aiPi + X aiPi + 2a r p r > a r 2p r - ^ Pi 

i>i* i>i*,ij^r i^r,i* y i>i* ,ij^r 

= a r (2p r - ( : Pi * -p r ) - (1 -Pr-Pi*)) = a r (4p r - 1) > 0. 

Thus we conclude that the second term of (7) is nonnegative and (5) holds. 



□ 


2.3 A approximation algorithm for monotone /c-submodular functions 

In this section, we show a polynomial-time randomized 2 k-i -approximation algorithm for maxi¬ 
mizing monotone /c-submodular functions. Our algorithm is described in Algorithm 3. We note 
that a similar algorithm and analysis appeared in [5] for the submodular welfare problem, which is 
a special case of the monotone /c-submodular function maximization problem. 

It is clear that Algorithm 3 runs in polynomial time. Below we consider the approximation 
ratio of Algorithm 3. 

Theorem 2.4. Let o be a maximizer of a monotone nonnegative k-submodular function f and let 
s be the output of Algorithm 3. Then E[/(s)] > 2 fc-i /(°)- 

Proof. By Lemma 2.2 it suffices to prove (4) for every 1 < j < n for c = 1 — p For simplicity of 
the description we shall omit the superscript (j) if it is clear from the context. 

We first consider the case (5 = 0. Since / is monotone, we have m = ai = 0 for all 1 < i < k. 
Hence, (4) clearly holds with c = 1 — p 
Now suppose f3 > 0. Our goal is to show 

x ( 1 -^) x y\ +1 - ( 9 ) 

l<i<k l<i<k 

If k = 1, then (9) follows since i* = 1 and both sides are equal to zero. Hence we assume k > 2. 
Let 7 = {k — 1) t = ft. Since / is a monotone /c-submodular function, we have that a* > 0 for all 
i € {1,..., k}. Then, we have 

X y^ ai * ~ a *) - X y i ai * - X yt i yi * = ~ y Vi * ' X yt i 

i^i* i^i* ^ i^i* 
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( 10 ) 





Algorithm 3 

Input: A monotone fc-submodular function / : {0,1,... , k} v —>• M+. 
Output: A vector s G {0,1,... , k} v . 

S i — 0. 


t <— k — 1 . 

for each e € V do 

Vi <- A £ii (s) for 1 < i < k. 

..t 


k 

£ v\ 

i =1 

if /3 / 0 then 


P 


Pi 

else 


(1 < i < k). 


Pi 


(1 <i< k ). 


I 1 if i = 1 
I 0 otherwise 
Let s(e) G {1 ,,k} be chosen randomly, with Pr[s(e) 


return s 


i] = pi for all * G {1,, k}. 


1 t 

From the weighted AM-GM inequality, a t + 1 * * b t + 1 < j^a + j^b holds for all a, b > 0. By setting 

a = ( 7 Ui*) t+1 and b = we have 


< 10 +(+<™*>‘ +1 + A 




t +1 
t \ t 


(ii) 


t+i t 


From Holder’s inequality, £, : ai < (£j«j* ) t+1 (£* l i+1 ) t+1 holds for any non-negative aj’s. By 
setting a,i = y\, we have 


( 11 ) <- 


1 ( 1 


7 \t + 1 


(79,-) ,+i +At(e»; +i )'(e 


mV 








7 


t + 1 


£<r-H)£ 


,t +1 


Thus we established (9) and we have k/(2k — l)-approximation by Lemma 2.2. 


□ 


3 Inapproximability 

As we remarked in the introduction, for a symmetric submodular function / : 2 V —> M + , a function 
g : {0, 1 ... , k} v —»• M + defined by 

k 

g{X u ...,X k ) = Y J f{Xi) ( {-V,..., X k } G { 0 ,... , k} v ) 

2=1 


9 












is /c-submodular. Hence one can consider an approximation algorithm for maximizing / by ap¬ 
plying an a-approximation algorithm for /c-submodular functions to g and then returning Xj £ 
argma x{f(Xj) \ j £ [A:]} for output (Xi,..., X k ) of the approximation algorithm. Let (X*,..., Xf) 
be a maximizer of g and X* be a maximizer of /. Since / is symmetric, we have g(X j",..., XT) > 
2 /(X*). Therefore we have kf(Xi) > Zj f(Z) = 9& i,- • ■,**)> ag(Xf ,.. .,X* k ) > 2 af(X*). 
Thus it gives a 2a//c-approximation algorithm for the symmetric submodular function maximiza¬ 
tion. 

It was proved by Feige, Mirrokni, and Vondrak [ 6 ] that any approximation algorithm for sym¬ 
metric submodular functions with polynomial queries cannot achieve the approximation ratio better 
that 1/2. This implies that the best approximation ratio for the /c-submodular maximization prob¬ 
lem is at most a < k/4. This argument, via embedding of a symmetric submodular function to 
a /c-submodular function, gives the tight approximation bound for bisubmodular function, but for 
k > 4 it does not give a nontrivial bound. 

Instead of embedding submodular functions to /c-submodular functions, in this section we shall 
directly extend the argument of [ 6 ] and establish the following bound. 

Theorem 3.1. For any e > 0, a + e )-approximation for the monotone k-submodular function 

maximization problem would require exponentially many queries. 

Proof. For simplicity we assume that £ is rational. Let V be a finite set with n = \V\ such that 
en is an integer. The framework of the proof is from [ 6 ] (see also [14]) and it proceeds as follows. 
We shall define a /c-subnrodular function / and a /c-submodular function gp for each /c-partition 

V = {A \,..., A k } of V, where a /c-partition means a partition of V into k subsets. Those functions 
look the same as long as queries are “balanced” (whose definition will be given below). Suppose 

V is randomly taken in the sense that each element is added to one of the k parts uniformly at 

random. Then it turns out that with high probability all queries would be balanced as long as 

the number of queries is polynomial in k and n. In particular, we cannot get any information 

about V. Thus one cannot distinguish / and gp by any deterministic algorithm with a polynomial 

number of queries. Hence, by Yao’s min-max principle, any (possibly, randomized) algorithm with 

a polynomial number of queries cannot distinguish / and gp and cannot achieve an approximation 

ratio better than maXiJ \ , which will be ^rf. 

max® gv{ x ) 7 ^ 

Now we define / and gp. For x = (X \,..., X k ) € {0,1,..., k} v , let uq{x) = |V \ Ui=i Xf\. We 
define f : {0,..., k} v —> Z + by 

f(x) = (k + 1 + 2 ke)n 2 — (k — l)no(ac ) 2 — 2(1 + ks)nno(x) (® € { 0 ,..., k} v ). 

To define gp , take any /c-partition V = {A \,..., A k }. For x = (X\, ..., X k ), let Cij(x) = \Xi n Aj\ 
for 1 < i < k and 1 < j < k, and let dj(x) = Zi=i °i,j+i -1 f° r 1 < j < fc, where the index is taken 
modulo k (0 is regarded as k). Then gp : {0,... , /c}' —> Z + is defined by 

9V (*) = /(*) + h v b ( x ) (x£{0,...,k} v ), 

l<a<b<k 

where 


hfp b (x) := (max{|d a (cc) — db(x)\ — en, 0}) 2 . 
The properties of / and gp are listed in the following claims. 
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Claim 3.2. f and gp for each k-partition V of V are nonnegative monotone k-submodular func¬ 
tions. 

Proof. Clearly they are nonnegative. To see the nronotonicity and fc-submodularity of /, take any 
x and e ^ supp(a;). Then 


A e ,if( x ) = 2(k — l)no(aj) + 2(1 + ke)n — 1. 

This is clearly nonnegative and hence / is monotone. Also, since no(-) is non-increasing, / is orthant 
submodular and hence / is /c-submodular by Theorem 1.1. 

Next we consider gp. Take any x and e ^ supp(cc), and suppose that e € Aj. When adding e 
into Xi , Cij(x ) increases by one and hence dj-i+ifx) increases by one. Hence for 1 < a < b < k 
we have 


A e ,ihp b {x) 


2((d a (x) - d b (x)) - en) + 1 
-2((d b (x) - d a (x )) - en) + 1 
< -2((d a (x) - d b (x)) - en) + 1 
2 ((d b (x) - d a (x)) - en) + 1 
0 


if a = j — i + 1 and d a (x) — d b (x) > en 
if a = j — i + 1 and d b (x) — d a (x) > en 
if b = j — i + 1 and d a (x) — d b (x) > en 
if b = j — i + 1 and d b (x) — d a (x) > en 
otherwise. 


Hence 


A £)i gp(x) = 2(k - l)n 0 (x)+ 2(1 + ke)n - 1 + ^ A e ^h! p l+l ' b (x). 

1 <b<k,b^j-i-\-l 


To see the orthant submodularity of gp, observe that 2no + A e^hip 4+1,6 is non-increasing for each 
b with b 7 ^ j — i + 1. Since 

A e ,igp(x) = 2(1 + fce)n - 1 + ^ ( 2 n 0 (a;) + A eii h^ l+1,b (x)), 

A e ^gp is non-increasing, implying the orthant submodularity. 

To see the monotonicity, let B + = {b \ dj-i + \(x) — d b (x) > en} and = {b \ d b (x) — 
dj-i+i(x) > en}. Note also that Ei < b <k d b( x ) = Ei< a <M<&<fc c a ,b( x ) = Ei< a <fc,i<6<fc \ X oPA\ < 
n. Hence 


^ e,i9V (®) 

= 2(k — l)n 0 (a;) + 2(1 + ke)n - 1 + ^ 2(dj_ i+ i(a;) - d b (x)) + (|H_| - \B + \)en + k - 1 

6SB+UB- 

> 2 I n — d b (x) I + 2 (k — |H + |)en > 0. 

\ 6SB+UB_ J 

This completes the proof of Claim 3.2. □ 

Claim 3.3. max^ f(x) = (k + 1 + 2 ke)n 2 and ma x x gp(x) > 2kn 2 (l — 0(e)). 


11 



Proof. Since / is fc-submodular, the maximum is attained for a ^-partition by Proposition 2.1, i.e., 
no(x ) = 0. Hence the maximum value of / is (k + 1 + 2 ke)n 2 . 

To see the second statement, take x = (X \...., X^) such that X t = Aj for 1 < i < k. Then 
d\(x) = n and dj(x ) = 0 for 2 < j < n— 1, and thus max^ g-p(x) > (k + l)n 2 + (A;— l)(n —en) 2 . □ 

Now take a random ^-partition V of V, and consider any deterministic algorithm that tries to 
distinguish / and g-p, where the algorithm do not know V. The algorithm issues some queries to 
the value oracle. Call a query to f(x) unbalanced if | di{x) — dj(x)\ > en for some i,j G [k], and 
otherwise balanced. Note that di(x) — dj[x) can be seen as a sum of independent random variables 
{Z e } ee y, where Z e = 1 if di(x) is increased due to e, that is, e € X^ n A i+ ^_i for some k, Z e = — 1 
if dj(x) is increased due to e, that is, e € X^ n for some k, and Z e = 0 otherwise. By 

ne 2 

Hoeffding’s inequality 1 , the probability that | di(x) — dj(x)\ > en for a query is at most 2e a - . 

2 

Hence by the union bound the probability that a query is unbalanced is at most k 2 e~~. Therefore, 

ne 2 

for any fixed sequence of e 4 queries, the probability that a query is unbalanced is still at most 
2 2 2 2 

k 2 e — = k 2 e Therefore, with probability at least 1 — k 2 e~~, all the queries will be 

balanced. As long as queries are balanced, the algorithm gets the same answer regardless of V, 
and it will never find out any information about the fc-partition V. In other words, with high 
probability, the algorithm will never distinguish between / and g-p. However, the maximum of / is 
(1 + k + 2 ke)n 2 while the maximum of gp is at least 2fc(l — 0(e))n 2 . This means that there is no 
polynomial-query algorithm with approximation ratio better than □ 

4 Approximation algorithms for skew-bisubmodular functions 

In this section, we discuss the problem of maximizing an a-bisubmodular function. An adaptation 
of the greedy algorithm is shown to achieve the approximation ratio of for a G [0,1]. 

This ratio converges to zero as a goes to zero. In order to improve the performance for small a, 
we give another simple approximation algorithm that achieves the approximation ratio of 
in Section 4.2. By taking the maximum of the outputs of these two algorithms, we obtain the 
approximation ratio of for any a G [0,1] (the minimum of the two ratio is achieved when 

a = jg). 

Concerning the maximum of an a-bisubmodular function, we have the following counterpart of 
Proposition 2.1. 

Lemma 4.1. For any a-bisubmodular function f : 3 1 —>• R_|_ with a G [0,1], there exists a partition 
of V that attains the maximum value of f. 

Proof. Suppose that ( S,T ) G 3 V attains the maximum value of /. By the a-bisubmodularity of /, 
we have 

af(S, V\S) + f(V\ T, T) > (1 + a)f(S, T), 

which implies that f(S,V\S ) = f(V \ T, T) = f(S, T). Thus the maximum value of / is attained 
by a partition of V. □ 

J Let Xi ,..., X n be independent random variables in [—1,1], and let X = 4(Xi + • ■ • + X n ). Then P(|X — E(A)| > 
t) < 2e 
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4.1 A randomized greedy algorithm 

We now extend the randomized greedy algorithm for the bisubmodular function (i.e., Algorithm 2 
for 2-submodular functions). Intuitively, a-bisubmodularity is a variant of bisubmodularity directed 
toward the first argument by parameter a. Following this intuition, we shall adjust the choice 
probability as shown in Algorithm 4. 


Algorithm 4 


Input: A nonnegative a-submodular function / : {0,1,2} V —>• M+. 
Output: A vector s G {0,1, 2} v . 
s<— 0. 

for each e G V do 
Vi A e ,if(s) for * = 1,2. 

if y 2 < 0 then 


Pi 


| 1 if * = 1 
[0 if * = 2. 
else if yi < 0 then 

f 0 if * = 1 
1 if * = 2. 


Pi 

else 

Pi 


ay 1 

ay 1 +y 2 

V2 


if i = 1 
if i = 2. 


I aj/i+j/2 

Let s(e) G {1,2} be chosen randomly, with Pr[s(e) 

return s 


i] = Pi for all i G {1, 2}. 


Note that, by the a-bisubmodularity of /, we have 


oA e> i/ (x) + A e , 2 f(x) > 0 


for any xG3 l and e ^ supp(a;), which implies ayi + yi > 0. 

The following theorem provides a performance analysis of this algorithm. 


( 12 ) 


Theorem 4.2. For any a G [0,1], the randomized greedy algorithm for maximizing a-bisubmodular 


functions provides an approximate solution within a factor of 


(l+v^) 2 ‘ 


Proof. Note that Lemma 2.2 holds for any function on {0,1,... , k}^ as long as the maximum 
is taken by a partition of V, which is the case in the a-bisubmodular function maximization by 
Lemma 4.1. Hence it suffices to prove (4) for every 1 < j < n for c = (l + a)/(2y / a). For simplicity 
of the description we shall omit the subscript ( j ). Our goal is to show 


(a*. - ai)pi + (a;* - a 2 )p 2 < c(yipi + y 2 p 2 ) 


(13) 


for each i* G {1,2} by using 


aai + a 2 > 0 

Vi > at (* = 1,2), 


13 








which follow from the a-bisubmodularity. 

If yi < 0, then a\ < 0 and y 2 > 02 > —aa± > 0, and hence cy 2 > y 2 > minjai —02, 0}, implying 
( 13 ). A symmetric argument also implies ( 13 ) if y 2 < 0. 

By ay\ + y 2 > 0, the remaining case is when 2/1 > 0 and y 2 >0. If i* = 1 , then 


(ai - a 2 )p 2 = 


< 


(ai - a 2 )y 2 < (1 + a)aiy 2 
<*2/i + 2/2 ~ ayi + 2/2 
(l + a)(a!2/i +2/1) = 1 + <* 
2 v / a(ci;2/i + 2/2) 2^/a 


< (1 + a)2/i2/2 
~ <*2/i + 2/2 

(2/1P1 + 2/2P2) ■ 


(1 + a)(ay 1 + y% - (y/ayi ~ 2/2) 2 ) 
c l y fa{ay\ + y 2 ) 


Thus ( 13 ) holds. On the other hand, if i* = 2 , 


{ „ 2 _ ol)p , = < A 




< 


<*2/i + 2/2 <*2/i + 2/2 <*2/i + 2/2 

(1 + <*)(<*2/i +2/1) l + <*, . x 

- o 7=7 -;-\ = 0 (yiPi + U2P2)- 

2y/a(ayi + y 2 ) 2 y/a 


(1 + a)(ay% + 2/2 - (\/<*2/i ~ 2/2) 2 ) 
2 v / a(ai/i + 2/2) 


Thus ( 13 ) holds. 


□ 


4.2 The second algorithm 

In this section, we describe another algorithm for maximizing a-bisubmodular functions, which 
achieves a better approximation ratio than the randomized greedy algorithm for small a. 

For an a-bisubmodular function / : 3 V —>• M+ with a € [0,1], we define f : 2 V -> R + 
by f'(X) = f(X, 0) for A € 2 1 . Since f{X) is a non-negative submodular function, we can 
apply the randomized double greedy algorithm of [2] to obtain a ^-approximate solution Z to the 
maximization of f. Our second algorithm for a-bisubmodular function maximization is rather 
simple: Take the better of (0, V) and (Z, 0). 

Theorem 4.3. The second algorithm for a-bisubmodular function maximization problem provides 
a - approximate solution for any a € [ 0 , 1 ]. 

Proof. Let (S, T) be an optimal solution. By the a-bisubmodularity, we have 

f(S, 0) + / (0, V) > f (0, 0) + af (0, T) + (1 - a)f(S, T), 
f (S, 0) + / (0, T) > f (0,0) + / (S, T), 

which imply 

(1 + a)f(S, 0) + / (0, V) > (1 + <*)/(0,0) + f(S, T ) > /(S, T). 

Let (A, B) be the output by the algorithm. Then, 

E[/(A B)] > max j 1/(S, 0), /(», V)} > |±A • )/(S, 0) + jA_/(0. V) 

TO «* + + '<*• a T >- D 

Combining this with the first algorithm, we obtain an approximate solution within a factor of 
max | , g -jL_. |. The minimum of this ratio is ^, which is achieved when a = 


14 
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